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ABSTRACT 

Photometric UBVI CCD photometry is presented for NGC 188 and Berkeley 17. 
Color-magnitude diagrams (CMDs) are constructed and reach well past the main- 
sequence turn-off for both clusters. Cluster ages are determined by means of isochrone 
fitting to the cluster CMDs. These fits are constrained to agree with spectroscopic 
metallicity and reddening estimates. Cluster ages are determined to be 7.0 ± 0.5 Gyr 
for NGC 188, and 10.0 ± 1.0 Gyr for Berkeley 17, where the errors refer to uncertainties 
in the relative age determinations. These ages are compared to the ages of relatively 
metal-rich inner halo/thick disk globular clusters and other old open clusters. Berke¬ 
ley 17 and NGC 6791 are the oldest open clusters with an age of 10 Gyr. They are 
2 Gyr younger than the thick disk globular clusters. These results confirm the status 
of Berkeley 17 as one of the oldest known open cluster in the Milky Way, and its age 
provides a lower limit to the age of the Galactic disk. 

Subject headings: Galaxy: disk - Galaxy: formation - open clusters and associations: 
general — open clusters and associations: individual (NGC 188, Be 17) 


1. Introduction 

Galaxy formation and evolution theory remains one of the great outstanding problems in 
contemporary astrophysics. Although considerable progress has been made in this field over the past 
decades, there are many unanswered questions regarding the formation of galaxies like the Milky 
Way exists. Observations within the Milky Way galaxy can be used to probe galactic evolution. In 
particular, determining the relative ages of the different stellar populations in the Milky Way - 
the halo, thick disk, thin disk, and bulge — by dating open and globular stellar clusters provides 
significant insight into the chronology of Galaxy formation (Liu & Chaboyer 2000; Salaris et al. 
2004). 

Owing to the great importance of the old open cluster population in probing the chemical 
and dynamic evolution of the Galaxy, considerable effort has been devoted to determining the 
physical parameters of these clusters (e.g. Phelps et al. 1994; Janes &: Phelps 1994; Scott et al. 
1995). However, relatively few of the oldest open clusters have received commensurate attention 
in astrometric, photometric, and spectroscopic studies, despite the great promise they hold in 
determining the age of the Galactic disk (notable exceptions include NGC 188 and NGC 6791). 
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This paper presents accurate photometric UBVI photometry for old open clusters NGC 188 
and Berkeley 17. A number of excellent studies have previously been performed for NGC 188 that 
have accurately determined its physical parameters (von Hippel & Sarajedini 1998; Sarajedini et 
al. 1999; Platais et al. 2003; Michaud et al. 2004; Stetson et al. 2004; VandenBerg Sz Stetson 2004). 
NGC 188 is used as a fiducial cluster, and a precise, relative age is obtained for Berkeley 17 relative 
to NGC 188. Additionally, by directly comparing it with previously published data, the NGC 188 
photometry serves as an independent test of the accuracy of the photometric calibration from the 
instrumental to the standard system. 

Since the first photometric study of NGC 188 by Sandage (1962) demonstrating that it belonged 
to the oldest open clusters in the Galactic disk, the cluster has been the subject of numerous studies. 
As the highest priority cluster in the WIYN Open Cluster Study, excellent, multi-color CMDs and 
proper-motion data have been published for this cluster (von Hippel & Sarajedini 1998; Sarajedini 
et al. 1999; Platais et al. 2003). In their UBVRI photometric study, Sarajedini et al. (1999) found 
that NGC 188 has an age of 7.0 ± 0.5 Gyr, a reddening of E(B-V) = 0.09 ± 0.02, and a distance 
modulus of (m-M)y = 11.44 ± 0.08. These values are in general agreement with other recent 
photometric studies of the cluster (Caputo et al. 1990; Twarog & Anthony-Twarog 1989). Stetson 
et al. (2004) obtained new data on NGC 188 and did a comprehensive review of existing data in 
the literature resulting in a large, homogeneous photometric database. These data were used by 
VandenBerg & Stetson (2004) to determine an age of 6.8 ± 0.7 Gyr for NGC 188. 

Berkeley 17 was discovered by Setteducati & Waever (1962), and is located at a low Galactic 
latitude in the direction of the Galactic anticenter. As a consequence of the cluster’s location as 
well as its large distance, the held in the direction of Berkeley 17 is highly reddened and greatly 
contaminated by held stars. Since the extensive open cluster study by Phelps et al. (1994), a 
general agreement has prevailed that Berkeley 17 is indeed the oldest known open cluster in the 
Galaxy. Soon after the Phelps et al. study, Kaluzny (1994) performed the first BVI photometry 
of the cluster, and established that Berkeley 17 is as old or somewhat older than the previously 
oldest known open cluster NGC 6791, based on comparisons of the morphologies of the two clusters’ 
CMDs. However, the Kaluzny study was hampered by poor weather, thereby preventing absolute 
determinations of Berkeley 17’s physical parameters. 

The hrst age determinations for Be 17 based on isochrone fitting were made by Phelps (1997) 
who found an age of 10-13 Gyr, a metallicity of -0.30 < [Fe/H] < 0.00, a reddening of 0.52 < 
E(B-V) < 0.68 and 0.61 < E(V-I) < 0.71, and a distance modulus of (rn-M)y = 14.05 ± 0.25. 
Although other attempts at establishing the age of Berkeley 17 have been carried out (Carraro et 
al. 1999b; Salaris et al. 2004), no further observational data has been published for the cluster. 

Section 2 describes the observations and the data reduction process, and §3 presents cluster 
CMDs. Section 4 describes our stellar models, isochrone fits, and the age determinations of the 
clusters. These ages are compared to the ages of other old open clusters and relatively metal-rich 
globular clusters in §5. Section 6 evaluates these results in the context of the formation of the 



- 3 - 


Milky Way. 


2. CCD Photometry 
2.1. Observations 

Observations were made at the MDM Observatory 1.3m McGraw-Hill telescope for thirteen 
nights between September 7 and 19, 2004. The Templeton camera, a thinned, backside-illuminated 
SITE 1024 x 1024 pixel CCD, was used, producing an image scale of 0.50 arcsec per pixel and a 
held of view of 8.5 x 8.5 arcmin. The observations are summarized in Table 1. 

Preliminary processing of the data was carried out within the IRAF 1 data reduction software 
and consisted of standard zero-level correction and hat helding. For the zero-level correction, 
between 15 and 20 bias frames were obtained throughout each night. Twilight and dawn hats were 
acquired every evening and used for hat helding. 

Four standard star helds from Landolt (1992) — PG1633, PG2213, and two helds in L95 - 
were observed throughout the potentially photometric nights. These helds consisted of a total of 
22 stars with a color range of —0.22 < (B — V) < 2.00 and were observed at airmasses ranging 
from 1.1 to 1.8 to ensure accurate determination of the photometric extinction and color terms. 


2.2. Standard Star Frame Reduction 

Aperture photometry was completed on the standard star frames using the APPHOT package 
within IRAF. The transformation equations between the instrumental and standard systems were 
established using the PHOTCAL package; these equations contain a zero point, an extinction term, 
and a linear color term, and are given by: 

u = U + u\ + u 2 X + U 2 ,{U — B) 
b = B + b 1 + b 2 X + b 3 (B -V) 
v = V + vi + v 2 X + v 3 (B - V) 
i = I + ii + i 2 X + 23(14 — I) 

where u, b , v. and i are the instrumental magnitudes, U, B , V, and I are the standard magnitudes, 
and X is the airmass. The color term coefficients represent a weighted mean average from the three 
photometric nights (September 10, 14, and 15), and were found to be 113 = —0.025, 63 = —0.040, 
v% = —0.014, and *3 = 0.014. The zero-point and extinction terms, however, showed sufficient 


1 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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night-to-night variation to be determined separately for each night, assuming the average color 
terms. No trends were observed in the residuals to the fits with respect to time, airmass, or color. 
Typical rms deviations of the fits to the standard values were 0.10 mag for U, and 0.02 mag for B, 
V, and I. 


2.3. Cluster Frame Reduction 

Point-spread function (PSF) photometry of the cluster frames was performed with the DAOPHOT 
II (Stetson 1992) photometry package. An average of 80 uncrowded stars from each frame were 
selected to construct a PSF. The PSF photometry was then carried out with the ALLSTAR pack¬ 
age by fitting the final PSF to detected profiles on each frame. The subtracted frames were then 
examined to ensure they contained only saturated stars, galaxies, and chip defects. Stars with 
magnitude errors larger than 0.1 were edited out of the final photometry file. 

Aperture corrections were determined by performing aperture photometry on 100-150 bright, 
uncrowded stars in each frame. The difference between the magnitudes given by the PSF and 
aperture photometry constituted the aperture correction, which was then applied to the original 
PSF photometry. Typical values for the aperture correction ranged between 0.1 and 0.2, and the 
values showed no trends in variation across the chip. 


3. Color-Magnitude Diagrams 

Color-magnitude diagrams (CMDs) for the clusters were constructed using the standard mag¬ 
nitudes given by the transformation equations (see Section 2.2). Stars detected in at least two 
frames of a given filter and exposure time were averaged, and were included in the final CMD if 
they also appeared in two other filters. The final Be 17 cluster CMDs are shown in Figures 1-2, 
and photometry files are given in Table 2 for NGC 188 and Table 3 for Be 17. Entries of 99.99 in 
these tables indicate that a star was not detected in that filter. 

The NGC 188 photometry includes only stars within the central region of the cluster, as the 
cluster has an estimated angular size of approximately 1°. The photometry extends around 6 
magnitudes fainter than its main-sequence turnoff at V ~ 15. The CMD includes 536 stars, and 
displays a well-defined main-sequence, subgiant branch, and red giant branch. A helium-burning 
clump cannot be accurately identified. A binary sequence located 0.75 magnitudes above the main 
sequence is detected, most notably in the ( V ., V — I) CMD. Our photometry covers a much smaller 
area, and is not quite as deep as the photometry presented by Sarajedini et al. (1999) and Stetson 
et al. (2004). Our photometry of NGC 188 was obtained so that a precise, relative comparison 
could be made be the well studied cluster NGC 188 and Be 17. In addition, the NGC 188 can be 
used to assess the absolute calibration of our photometry, via a comparison to other data sets. 
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This paper’s NGC 188 photometry was compared to that of Stetson et al. (2004) on a star 
by star basis. We were able to compare photometry of 533 stars in B, 531 stars in V, 525 stars 
in I and 304 stars in U. Zero-point offsets of +0.026 mag in U, +0.013 mag in B, +0.031 mag in 
V and +0.044 mag in I were found, in the sense that our photometry was fainter in all filters. 
The residuals between our photometry and Stetson et al. (2004) had no systematic trends with 
magnitude or colour. A sample comparison is shown in Figure 3, where the differences between 
our V photometry and Stetson et al. (2004) are plotted as a function of magnitude and colour. A 
similar comparison with the NGC 188 photometry of Sarajedini et al. (1999) determined somewhat 
smaller zero point offsets of +0.025 mag in U, —0.028 mag in B, —0.002 mag in V and +0.016mag, 
in the sense that our photometry was fainter in U and I, and brighter in B and V. Again, the 
comparison revealed no systematic trends with magnitude or colour. It is not clear whose data 
is best calibrated onto the standard system. Based upon these comparisons it appears that our 
photometry is on the Landolt (1992) photometric system to an accuracy of ±0.03 mag. 

As a result of Berkeley 17’s large distance and low Galactic latitude, its CMD is highly con¬ 
taminated by field stars. The CMD consists of a total of 1239 stars, and extends approximately 
4 magnitudes fainter than the main-sequence turnoff, located at V ~ 18. A highly contaminated 
main-sequence can be discerned, along with a well-populated subgiant branch and a scarcely pop¬ 
ulated giant branch. Neither a helium-burning clump not a binary sequence can be discerned. The 
(U — B,V — I) colour-colour diagram can be used to estimate reddening, and thereby remove many 
foreground stars (which typically have a lower reddening) from the CMD of Berkeley 17. The 
colour-colour of NGC 188 (which has a low reddening, and few foreground stars) was compared 
with the Be 17 colour-colour diagram to partially eliminate field star contamination (Figure 4). 
The cleaned Berkeley 17 CMD contains 1083 stars, and is shown in Figure 5. Stars which are 
deemed to be non-members of Bel7 are denoted by a zero in the final column of Table 3. 

A comparison between this paper’s (V, B — V) photometry (before field star removal) and that 
of Phelps (1997) is shown in Figure 6. While the Phelps photometry (which surveyed a somewhat 
smaller region) reaches approximately a magnitude deeper, this paper’s photometry displays both 
a tighter main-sequence and a more well-defined main-sequence turnoff. The Phelps photometry 
is redder than our photometry, with 5(B — V) ~ 0.1 mag. Given the good agreement of our NGC 
188 photometry with the photometry from Sarajedini et al. (1999); Stetson et al. (2004) we believe 
that our photometry is on the standard Landolt (1992) system to within ±0.03 mag. 


4. Isochrone Fitting and Age Determinations 
4.1. Stellar Models and Isochrones 

Stellar evolution tracks were constructed using our stellar evolution code (Chaboyer et al. 1999, 
2001) for masses in the range of 0.5 M 0 to 2.0 M 0 in increments of 0.05 M 0 . The models were 
evolved in 2000 time steps from the zero age main-sequence through the red giant branch. The 
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models include the diffusion of helium and heavy elements. A solar calibrated model was calculated, 
yielding Zq = 0.02 and Y & = 0.275 (initial abundances). For other metallicities, dY/dZ = 1.5 was 
assumed (corresponding to a primordial helium abundance of Ybbn = 0.245). The solar calibrated 
mixing length was used for all the models. The isochrones were created for metallicities of [Fe/H] 
= +0.20, +0.10, +0.00, -0.12, -0.30, -0.40, -0.60, -0.70, -0.82, and -1.00. 

The transformation from the theoretical log(L) versus log(T e yj) plane to observed colors and 
magnitudes for the constructed isochrones was completed using two separate color transformations, 
as this represents one of the greatest uncertainties of the theoretical isochrones. The first transfor¬ 
mation utilized the Kurucz (1993) model atmospheres as described in Chaboyer et al. (1999), and 
the second, the VandenBerg & Clem (2003) color tables. The isochrones were constructed in the age 
range of 100 Myr to 15 Gyr for each of the aforementioned metallicities and color transformations, 
thereby allowing for considerable flexibility in performing the isochrone fits to the cluster CMDs. 


4.2. Isochrone Fitting 

One of the most important isochrone input parameters is metallicity, and therefore, accurate 
cluster metallicities are crucial to obtaining reliable age estimates. We are interested in obtaining 
precise relative ages for NGC 188 and Be 17. For this reason, we elected to use the metallicity 
estimates from Friel et al. (2002), who conducted an extensive, homogeneous spectroscopic survey 
of a large number of old open clusters, seeking to accurately determine both cluster metallicities 
and reddenings. This paper’s clusters, NGC 188 and Berkeley 17, were both included in their inves¬ 
tigation. As the Friel et al. (2002) metallicities for these two clusters were done in a homogeneous 
manner, it is likely that the metallicity difference between the two clusters was robustly measured. 

For relative age determinations, using precise relative metallicity determinations is more im¬ 
portant than the zero-point of the metallicity scale. In this regard, high-resolution spectroscopy of 
M67 (Hobbs & Thorburn 1991; Tautvaisiene et al. 2000) found [Fe/H] values which are about 0.1 
dex higher than Friel et al. (2002). Thus, the Friel et al. (2002) study may be underestimating the 
metallicity of the clusters. This implies that our absolute age determinations could be systemat¬ 
ically in error, but does not effect our relative age determinations. A 0.1 dex increase in [Fe/H] 
would systematically decrease our age estimates by approximately 8%. 

Friel et al. (2002) used Fe I and Fe-peak blends to measure metallicities and the H [3 strength of 
main sequence stars spanning the distance to the cluster were used to determine reddening values. 
Friel et al. (2002) determined [Fe/H] = —0.10 ± 0.09 and E(B — V) = 0.04 ± 0.04 for NGC 188. 
However, Friel et al. (2002) note that their H (3 reddening value is lower than published estimates 
and adopt E(B — V) = 0.08. For Be 17, Friel et al. (2002) found [Fe/H] = —0.33 ± 0.12 and 
E(B-V) = 0.58 + 0.09. 

In fitting the isochrones to the CMDs of the two clusters, fits are made to the upper main- 
sequence, main-sequence turnoff and sub-giant regions. The colour of the red giant branch was 
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not used in the fitting process, as it’s location in the theoretical models depends sensitively on 
super-adiabatic convection which is not well understood. Similarly, the colour of the lower main 
sequence was not used in the fit, as there are difficulties with the conversion between theoretical 
temperatures and colours for cool main sequence stars. The range of acceptable isochrone fits is 
constrained by the requirement of simultaneous fits in the ( V., B — V) and (V, V — I) CMDs. The 
metallicity, reddening, and distance modulus are required to be identical in both CMDs, assuming 
E(V — I) = 1.25 x [E(B — V )] (Dean et al. 1978). Metallicities and reddenings are required to be 
consistent with those obtained by Friel et al. (2002). Isochrone fit parameters are summarized in 
Table 4. 


4.3. Discussion 

4.3.1. NGC188 

As a result of the high quality of the photometry, and its distinct CMD morphology NGC 
188 is determined to be 7.0 ± 0.5 Gyr old. The results obtained using the two individual sets 
of color transformations are consistent, and an excellent agreement with the spectroscopic values 
from Friel et al. (2002) is also evident. A sample isochrone fit is shown in Figure 7. An age of 
7.0 ± 0.5 Gyr agrees well with the previous age determination of 7.0 ± 0.5 Gyr, by Sarajedini et 
al. (1999); 6.8 ± 0.7Gyr by VandenBerg & Stetson (2004) and 6.4Gyr by Michaud et al. (2004). 
These studies utilized different photometry and stellar evolution models. The agreement between 
the different studies suggest that the age of NGC 188 has been robustly determined. 


4-3.2. Berkeley 17 

The age of Berkeley 17 is found to be 10.0 ±1.0 Gyr. As a result of the larger reddening 
uncertainty, greater field star contamination, and the less distinct cluster CMD morphology, the 
established age for Berkeley 17 has a greater uncertainty than does NGC 188. However, by using 
the constraints on reddening and metallicity from spectroscopy and requiring simultaneous fitting 
in (V, B — V) and (V, V — I), reliable results are obtained. Similar ages are obtained using the two 
separate color transformations. A sample isochrone fit is shown in Figure 8. 

Our age for Be 17 is on the low side of the estimates made by Phelps (1997), who determined 
an age of 12.5 Gyr based on the VandenBerg (1985) isochrones, and 10.0-13.0 Gyr based on the 
Bertelli et al. (1994) isochrones. However, as no accurate spectroscopic study had been performed 
for Berkeley 17 at the time, Phelps (1997) allowed for a large range of metallicities and reddenings. 
Additionally, simultaneous fitting in the (V, B — V), ( V ., V — I) CMDs was not an important factor 
in Phelps’ isochrone fitting process, mainly because the isochrones were not normalized to the Sun 
either in terms of helium abundance or in terms of the mixing-length parameter. 



On the other hand, Carraro et al. (1999a) found an lower estimate of 9.0 =h 1.0 Gyr using the 
Girardi et al. (2000) isochrones on both the Kaluzny (1994) and Phelps (1997) photometry. Carraro 
et al. (1999b), using the same methodology as Carraro et al. (1999a) made age determinations of 
previously well-studied clusters NGC 188 and NGC 6791, for which age estimates of 6.5 ± 0.5 Gyr 
(NGC 188) and 8.5 ± 0.5 Gyr (NGC 6791) were determined. The Carraro et al. (1999b) results 
imply an age difference between NGC 188 and Berkeley 17 of approximately 2.5 Gyr, consistent 
with the results obtained in this paper. The fact that our age estimates for Be 17 agree with Phelps 
(1997) and Carraro et al. (1999a) is surprising, given that our photometry has a significant offset 
from the Phelps (1997) photometry (see Figure 6). 


5. Comparison to Other Clusters 

In order to study the early formation history of the Milky Way, it is useful to compare the 
derived ages of NGC 188 and Bel7 to other old open clusters and globular clusters. As accurate es¬ 
timates of the metallicity of a cluster are critical for age determinations, open clusters were selected 
from the Friel et al. (2002) metallicity study which included [Fe/H] measurements for 39 clusters. 
From this list, clusters were selected which had good quality photometry and a morphological age 
index (Janes & Phelps 1994) greater than 2.0 (corresponding to an age of approximately 1 Gyr). 
The final sample of open clusters contains 18 clusters. 

Globular clusters were selected from the extensive, homogeneous HST photometric study of 
globular clusters carried out by Piotto et al. (2002). We elected to restrict our sample to relatively 
metal-rich ([Fe/H] > —1.2) clusters located in the inner part of the galaxy ( Rgc < 8kpc). Requiring 
the globular clusters to be relatively metal-rich makes the age comparisons to the more metal-rich 
open clusters more direct. The restriction on Galacto-center distance was made in an attempt to 
only include clusters which were formed as part of the initial Galactic collapse, and not as part of a 
latter accretion event. Salaris & Weiss (2002) found that globular clusters within a Galactocentric 
radius of 8 kpc are nearly coeval, implying that the inner halo globular clusters formed in the initial 
Galactic collapse, whereas the outer halo — where clusters were found to have a larger range in 
age — was likely formed through the accretion of extra-galactic fragments. Finally, the globular 
classification provided in the recent study by Mackey & van den Bergh (2005) was used to exclude 
clusters belonging to the young halo (YH), as these are believed to have been formed in external 
satellites eventually accreted by the Milky Way. Our globular cluster sample contains 13 clusters 
with B,V photometry from Piotto et al. (2002). For three of these clusters — 47 Tuc, NGC 5927, 
and NGC 6652 — HST VI photometry was also available. 

Given the sensitivity of isochrone fitting to the assumed metallicity, it is important that a 
consistent metallicity scale be used for the open and globular clusters. While numerous studies 
have attempted to determine uniform globular cluster (Carretta &; Gratton 1997; Rutledge et al. 
1997; Kraft & Ivans 2003) and open cluster (Piatti et al. 1995; Twarog et al. 1997; Friel et al. 2002) 
metallicities, these studies have generally not overlapped between the two cluster populations. 
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Rutledge et al. (1997) used equivalent width observations to establish a linear relationship 
between Ca II triplet equivalent widths and the globular cluster abundance scale of Carretta & 
Gratton (1997). This lead to a relative metallicity abundance ranking for a total of 71 globular 
clusters. Cole et al. (2004) sought to extend the Ca II triplet abundance scale to the younger ages 
and higher metallicities of open clusters. New spectra were obtained for a sample of five globular 
clusters and six old open clusters, and using globular cluster abundances on the Carretta & Gratton 
(1997) scale and open cluster abundances on the Friel et al. (2002) scale, a linear relationship 
between [Fe/H] and Ca II line strength was found. A crucial point in their study was to establish 
the compatibility of the two metallicity scales, given that the calibration of two scales differs in 
two respects: the atmosphere models utilized, and the effective temperature derivations. Cole et 
al. (2004) conclude that the systematic offsets between the two scales are comparable to or smaller 
than their respective internal uncertainties. They believe that the offset resulting from the choice 
of model atmospheres is of similar magnitude but opposite sign of that of the choice of effective 
temperature scale. Cole et al. (2004) conclude that the offset between the two metallicity scales is 
no more than 0.1 dex. Globular and open cluster metallicities were thus taken from Rutledge et 
al. (1997) and Friel et al. (2002). 

Individual ages for all of the clusters were determined using main sequence fitting, in the same 
manner that was done for NGC 188 and Pal 11. In making these fits, the metallicity of the isochrone 
was required to be consistent to within 0.1 dex of the cluster metallicity. The isochrones were 
constructed using scaled solar compositions, while globular clusters are typically enhanced in their 
a-element (O, Mg, Si, S, and Ca) abundances. To take this into account, the relationship between 
the heavy elements mass fraction ( Z ) in the isochrones and [Fe/H] followed the relationship given 
by Chieffi et al. (1991) and Chaboyer et al. (1992). In making this correction, globular cluster 
[a/Fe] abundances were assumed to follow the relationship between [a/Fe] and [Fe/H] found in 
Edvardsson et al. (1993). The a-abundances of 47 Tuc were found to correspond well with this 
relationship (Carretta et al. 2004). The derived ages are listed in Table 5. 

The globular clusters in the sample can be divided up into three categories based on their 
location in the Galaxy and metallicity: thick disk, bulge, and inner halo. The three thick disk 
globular clusters, 47 Tuc, NGC 5927, and NGC 6838 are found to have a mean age of 12.3 ± 0.4 
Gyr. The three bulge clusters, NGC 6304, NGC 6624, and NGC 6637 have a mean age of 12.7±0.2 
Gyr. The seven clusters in the inner halo have an age range of 11.8-14.0 Gyr, with a mean age of 
13.0 ±0.9 Gyr. 2 Collectively, the thirteen high-metallicity globular clusters in the sample are found 
to be effectively coeval, with an average age of 12.8 ± 0.7 Gyr and with all ages falling within 2a 
of the mean. We note that in general, our globular cluster ages are somewhat older (~ 10%) than 
other studies which attempt to get good absolute ages for globular clusters (Salaris & Weiss 2002; 


2 The errors in the mean ages of the three globular cluster populations represent the standard deviation of the 
ages with respect to the mean, and do not include errors due to isochrone fitting. The errors in the ages of individual 
clusters given in Table 5 represent relative, rather than absolute, errors, and thereby, no conflict with cosmological 
constraints on the age of the Universe exists. 
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Krauss & Chaboyer 2003; Lewis et al. 2005). 

The oldest open clusters are found to be NGC 6791, Berkeley 17, and Collinder 261, with ages 
of 10.0 ± 1.0 Gyr, 10.0 ± 1.0 Gyr, and 8.0 ± 1.0 Gyr, respectively. As a result of its extreme physical 
parameters — it is one of the oldest, most metal-rich, and most populous open clusters — NGC 
6791 has been investigated in countless photometric and spectroscopic studies. The age obtained 
for NGC 6791 in this study, 10.0 ±1.0 Gyr, is in excellent agreement with the result found by 
Salaris et al. (2004), who determined an age of 10.2 ± 1.2 Gyr. 

Salaris et al. (2004) used a combination of main sequence fitting and a new calibration of 
the morphological age indicator to determine the age of 71 old open clusters. The age of NGC 
188 from isochrone fitting was found to be 6.3 ± 0.8 Gyr, while the derived age for Be 17 was 
10.1 ± 2.8 Gyr. Our new photometry has allowed us to put a much tighter constraint on the age of 
Be 17. A comparison between the 17 open cluster ages in common between this study and Salaris 
et al. (2004) is found in Figure 9. In general, our open cluster ages are found to be consistent, 
within the error bars, with the ages obtained in Salaris et al. (2004). On average, the error bars 
in our age determinations are ~ 50% smaller than the error bars in the Salaris et al. (2004) age 
determinations. 

The only cluster for which our age disagrees with the age obtained by Salaris et al. (2004) is 
Berkeley 20. Be 20 is one of the most distant known open clusters, with a Galactocentric distance of 
16.4 kpc. The BVI photometry for this cluster Durgapal et al. (2001) contains significant field star 
contamination. The main-sequence turnoff, the most important region in the CMD in performing 
the isochrone fitting, was poorly-defined. However, it is evident from the cluster CMD that no so- 
called “hook” exists in the main-sequence turnoff region, a characteristic CMD feature of clusters 
with ages below around 5.0 Gyr. Thereby, it is unlikely the cluster has an age below 5.0 Gyr, as 
suggested by the Salaris et al. (2004) study. Furthermore, simultaneous fitting in the ( V ., B — V) 
and (V, V — I) CMDs was only possible with isochrones of around 6 Gyr, for which the reddening 
value was also consistent with the value given in the literature ( E(B — V) = 0.14). Berkeley 20 is 
thus determined to be 6.0 ± 1.0 Gyr old, which is consistent with the age of 5.0 Gyr obtained from 
isochrone fitting by Durgapal et al. (2001) 


6. Summary 

This paper confirms the well-established age of old open cluster NGC 188 at 7.0 ± 0.5 Gyr. Be 
17 is found to be 10.0 ± 1.0 Gry old. As a result of the internal consistency of the ages obtained 
for the two clusters, the age difference of Be 17 relative to NGC 188 is highly robust: Be 17 is 
3.0 ± 1.1 Gyr older than NGC 188. These ages were determined using the metallicity scale of Friel 
et al. (2002). There are some indications that the Friel et al. (2002) [Fe/H] measurements are too 
low by about 0.1 dex as high resolution spectroscopy studies of M67 (Hobbs & Thorburn 1991; 
Tautvaisiene et al. 2000) find [Fe/H] values 0.1 dex higher than Friel et al. (2002). If the Friel et al. 
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(2002) metallicities are 0.1 dex too low, then our ages need to be revised downward by about 8%. 

Using the same age determination technique, the Friel et al. (2002) metallicity measurements 
and photometry from the literature, ages of 16 other old, open clusters were determined. Be 17 
was found to be the oldest open cluster, with the same age as the well studied cluster NGC 6791. 
The ages of 13 relatively metal-rich globular clusters were determined using the same methodology 
as the open cluster ages. Table 5 summarizes the age determinations. These ages were determined 
using isochrone fitting in a self-consistent manner. This leads to precise relative ages, but the 
absolute ages have considerably larger uncertainties. We note that the globular cluster ages we 
derive are ~ 8% larger than the accurate absolute ages found by Krauss & Chaboyer (2003). This 
is the same age reduction estimated for the open clusters, based upon the suggestion that the Friel 
et al. (2002) metallicity values are too low by 0.1 dex. 

The globular clusters have ages in the range 11.8 - 14.0 Gyr, with no evidence for an age range. 
The thick disk globular clusters (the most metal-rich clusters in our sample) were found to have 
an average age of 12.3 ± 0.4 Gyr, while the two oldest open clusters in our sample have an age of 
10.0 ±0.7 Gyr. Thus, the oldest open clusters in the thin disk are found to be 2.3 ±0.8 Gyr younger 
than the thick disk globular clusters. 

In contrast, Salaris et al. (2004) find no significant age difference between the thin and thick 
disks in the Galaxy. However, their age determinations generally had larger error bars than our 
determinations, and Salaris et al. (2004) determined the age of two thick disk globular clusters, 
while we determine the age of three thick disk globular clusters. Combined, these two effects made it 
difficult for Salaris et al. (2004) to determine if age differences less than 3 Gyr existed. In summary, 
we find that the oldest open clusters imply that the thin disk started forming 2.3 ± 0.8 Gyr after 
the formation of the thick disk. 


We would like to thank the referee, whose thoughtful reading of our manuscript lead to several 
improvements. Research supported in part by a NSF CAREER grant 0094231 to BC. BC is a 
Cottrell Scholar of the Research Corporation. 
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Table 1. Log of Observations 


Date 

Cluster 

Exposure (s) 

Filter 

Airmass 

FWHM (arcsec) 

2004 Sep 10 

Be 17 

3 

X 

400 

U 

1.4 

1.8 



3 

X 

90 

B 

1.3 

1.7 



3 

X 

40 

V 

1.4 

1.4 



3 

X 

40 

I 

1.4 

1.4 


Be 17 Off-Field 

3 

X 

400 

u 

1.2 

1.8 



3 

X 

90 

B 

1.1 

1.8 



3 

X 

40 

V 

1.2 

1.4 



3 

X 

40 

I 

1.2 

1.4 

2004 Sep 15 

NGC 188 

3 

X 

120 

u 

1.7 

2.5 



2 

X 

1200 

u 

1.7 

2.5 



3 

X 

90 

B 

1.7 

2.0 



2 

X 

900 

B 

1.7 

2.0 



3 

X 

60 

V 

1.7 

2.0 



2 

X 

600 

V 

1.7 

2.0 



3 

X 

30 

I 

1.7 

2.0 



2 

X 

300 

I 

1.7 

2.0 


Be 17 

2 

X 

1200 

u 

1.2 

1.8 



2 

X 

900 

B 

1.1 

1.8 



2 

X 

600 

V 

1.3 

1.6 



2 

X 

300 

I 

1.3 

1.7 

2004 Sep 16 

Be 17 Off-Field 

2 

X 

1200 

u 

1.1 

2.0 



2 

X 

900 

B 

1.1 

1.7 



2 

X 

600 

V 

1.4 

1.6 



2 

X 

300 

I 

1.2 

1.5 
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Table 2. NGC 188 Photometry File Sample 


Number 

X (pix) 

Y (pix) 

U (mag) 

B (mag) 

V (mag) 

I (mag) 

i 

921.222 

1.588 

99.99 

19.92 

18.74 

17.31 

2 

265.934 

4.543 

99.99 

21.26 

19.83 

18.05 

3 

62.563 

5.192 

17.45 

17.18 

16.39 

15.57 

4 

12.556 

9.095 

99.99 

20.52 

19.17 

17.60 

5 

63.443 

9.635 

19.66 

18.52 

17.28 

15.94 

6 

405.515 

9.803 

17.64 

17.53 

16.83 

16.00 

7 

350.567 

11.787 

15.38 

15.39 

14.78 

14.03 

8 

314.874 

15.181 

19.74 

18.79 

17.72 

16.44 

9 

634.737 

17.929 

16.65 

16.47 

15.81 

15.04 
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Table 3. Berkeley 17 Photometry File Sample 


Number 

X (pix) 

Y (pix) 

U (mag) 

B (mag) 

V (mag) 

I (mag) 

member? 

i 

479.796 

1.618 

99.99 

21.98 

20.40 

17.79 

i 

2 

445.481 

2.832 

99.99 

20.69 

19.06 

16.80 

i 

3 

677.789 

3.111 

99.99 

21.37 

19.97 

18.11 

i 

4 

41.351 

5.832 

99.99 

21.27 

20.14 

18.57 

i 

5 

391.842 

5.954 

20.00 

19.13 

17.79 

99.99 

i 

6 

692.070 

5.973 

99.99 

20.86 

19.59 

18.04 

i 

7 

608.705 

6.829 

99.99 

21.66 

20.43 

18.79 

i 

8 

689.438 

7.807 

20.17 

19.92 

18.94 

17.57 

i 

9 

56.520 

9.221 

21.39 

20.78 

19.61 

18.00 

i 






Table 4. 


Isochrone Fit Parameters 


Cluster 

[Fe/H] 

E(B-V) 

(rn-M)y 

Age (Gyr) 

Color Transformation 

NGC 188 

+0.00 

0.06 

11.40 

6.5 

VandenBerg & Clem 


-0.12 

0.07 

11.35 

7.0 

VandenBerg & Clem 


-0.12 

0.07 

11.35 

7.0 

Kurucz 


-0.12 

0.08 

11.45 

6.5 

Kurucz 

Be 17 

-0.30 

0.59 

14.20 

9 

VandenBerg & Clem 


-0.30 

0.57 

14.10 

10 

VandenBerg & Clem 


-0.40 

0.61 

14.20 

9 

VandenBerg & Clem 


-0.40 

0.60 

14.10 

10 

VandenBerg & Clem 


-0.40 

0.58 

14.10 

10 

Kurucz 


-0.40 

0.56 

14.00 

11 

Kurucz 






Table 5. Cluster Parameters 


Name 

E(B - V) 

D 

Rgc 

Class 3, 

[Fe/H] 

Age (Gyr) 

Reference 

47 Tuc 

0.04 

4.5 

7.4 

TD 

-0.78 

12.8 

± 

1.0 

i 

NGC 5904 

0.03 

7.5 

6.2 

IH 

1.12 

12.0 

± 

1.0 

2 

NGC 5927 

0.45 

7.6 

4.5 

TD 

-0.64 

12.0 

± 

1.0 

2,3 

NGC 6171 

0.33 

6.4 

3.3 

IH 

-0.95 

13.5 

± 

1.5 

2 

NGC 6218 

0.19 

4.9 

4.5 

IH 

1.14 

14.0 

± 

1.0 

2 

NGC 6304 

0.53 

6.0 

2.2 

B 

-0.66 

13.0 

± 

1.0 

2 

NGC 6362 

0.09 

7.6 

5.1 

IH 

-0.99 

14.0 

± 

1.0 

2 

NGC 6624 

0.28 

7.9 

1.2 

B 

-0.70 

12.5 

± 

0.5 

2 

NGC 6637 

0.16 

9.1 

1.9 

B 

-0.78 

12.5 

± 

0.5 

2 

NGC 6652 

0.09 

10.1 

2.8 

IH 

-0.85 

11.8 

± 

1.0 

2,4 

NGC 6712 

0.45 

6.9 

3.5 

IH 

-0.94 

12.5 

± 

0.5 

2 

NGC 6723 

0.05 

8.7 

2.6 

IH 

-0.96 

13.5 

± 

1.5 

2 

NGC 6838 

0.25 

4.0 

6.7 

TD 

-0.73 

12.0 

± 

1.0 

2 

NGC 188 

0.08 

1.5 

9.3 

OC 

-0.10 

7.0 

± 

0.5 

5 

NGC 1193 

0.12 

4.0 

12.0 

OC 

-0.51 

4.0 

± 

0.5 

6 

NGC 2141 

0.33 

4.1 

12.4 

OC 

-0.33 

3.0 

± 

0.5 

7 

NGC 2204 

0.08 

4.0 

11.5 

OC 

-0.32 

2.0 

± 

0.5 

8 

NGC 2243 

0.06 

3.7 

10.8 

OC 

-0.49 

3.8 

± 

0.8 

9 

NGC 2420 

0.05 

2.3 

10.6 

OC 

-0.38 

2.0 

± 

0.5 

10 

NGC 2506 

0.07 

3.3 

10.9 

OC 

-0.44 

2.0 

± 

0.5 

11 

NGC 2682 

0.05 

0.8 

9.05 

OC 

-0.15 

5.0 

± 

0.5 

12 

NGC 6791 

0.15 

4.2 

8.12 

OC 

+0.11 

10.0 

± 

1.0 

13 

NGC 6819 

0.16 

2.3 

8.18 

OC 

-0.11 

2.5 

± 

0.5 

14 

NGC 6939 

0.50 

1.2 

8.70 

OC 

-0.19 

1.6 

± 

0.4 

15 

Be 17 

0.58 

2.7 

11.2 

OC 

-0.33 

10.0 

± 

1.0 

5 

Be 20 

0.14 

8.4 

16.4 

OC 

-0.61 

6.0 

± 

1.0 

16 

Be 21 

0.76 

5.0 

13.5 

OC 

-0.62 

2.0 

± 

0.5 

17 

Be 31 

0.20 

3.8 

12.0 

OC 

-0.40 

3.5 

± 

0.5 

18 

Be 32 

0.15 

3.1 

11.4 

OC 

-0.50 

5.0 

± 

1.0 

19 

Be 39 

0.12 

3.8 

11.5 

OC 

-0.26 

6.5 

± 

1.5 

8 

Cr 261 

0.33 

2.4 

7.5 

OC 

-0.16 

8.0 

± 

1.0 

20 


a Classification: B Bulge globular cluster; TD Thick Disk globular cluster; IH Inner Halo 
globular cluster; OC Open Cluster 


References. — Galactic Coordinates, reddenings and distances for globular clusters are 
from Harris (1996); Galactic Coordinates, reddenings and distances for open clusters are 
from Friel et al. (2002); Sources for photometry: (1) Kaluzny et al. (1998); (2) Piotto 
et al. (2002); (3) Feltzing & Gilmore (2000); (4) Chaboyer et al. (2000); (5) This paper; 
(6) Kaluzny (1988); ((7) Carraro et al. (2001); (8) Kassis et al. (1997); (9) Bergbusch et 
(al. (1991); (10) Anthony-Twarog et al. (1990); (11) Kim et (al. (2001); (12) Sandquist 
(2004); (13) Stetson et al. (2003); (14) (Kalirai et al. (2001); (15) Andreuzzi et al. (2004); 
(16) Durgapal et (al. (2001); (17) Tosi et al. (1998); (18) Guetter (1993); (19) (Kaluzny 
& Mazur (1991); 20) Gozzoli et al. (1996); All OC photometry (obtained from WEBDA, 
maintained by Jean-Claude Mermilliod (http://obswww.unige.ch/webda/navigation.html) 
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Fig. 4.— NGC 188 (left) and Berkeley 17 (right) (U — B,V — I) Color-Color Diagrams. Eliminated 
stars in the Berkeley 17 Color-Color Diagram are marked with open squares. 


U-B 
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Fig. 7.— Best simultaneous isochrone fits for NGC 188. Isochrones constructed with the Van- 
denBerg & Clem (2003) colour transformation with [Fe/H] = —0.12 for 6.5, 7.0, and 7.5 Gyr are 
plotted assuming E(B — V) = 0.07 and (m — M)y = 11.35. 
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Fig. 8.— Best simultaneous isochrone fits for Berkeley 17. Isochrones constructed with the Van- 
denBerg & Clem (2003) colour transformation with [Fe/H] = —0.30 for 9, 10 and 11 Gyr are plotted 
assuming E(B — V) = 0.57 and (m — M)y = 14.10. 
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Fig. 9.— Age determinations for the 17 clusters in common between this study (horizontal axis) 
and the ages derived by Salaris et al. (2004). The dashed line has unit slope and illustrates perfect 
agreement; it is not a fit to the data. The only cluster whose age, within the error bars, does not 
agree is Be 20 (with an age of 6.0 Gyr in our study). 




